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Abstract 



We present a measurement of the tt production cross section in pp collisions at yfs = 1.96 TeV 
using events containing a high transverse momentum electron or muon, three or more jets, and 
missing transverse energy. Events consistent with tt decay are found by identifying jets containing 
candidate heavy-flavor semileptonic decays to muons. The measurement uses a CDF Run II data 
sample corresponding to 2 fb -1 of integrated luminosity. Based on 248 candidate events with three 
or more jets and an expected background of 79.5 ± 5.3 events, we measure a production cross 
section of 9.1 ± 1.6 pb. 

PACS numbers: 13.85M, 13.85Qk, 14.65Ha 
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I. INTRODUCTION 



Top quark pair production in hadronic collisions in the standard model proceeds via ei- 
ther quark- ant iquark annihilation or through gluon-gluon fusion. At the Fermilab Tevatron 
collider, with a center-of-mass energy of 1.96 TeV, the production is expected to be domi- 
nated by quark- ant iquark annihilation. For a top mass of 175 GeV/c 2 the theoretical cross 
section is calculated to be 6.6 ± 0.6 pb [1] and decreases by approximately 0.2 pb for each 
1 GeV/c 2 increase in the top mass over the range 170 GeV/c 2 <M top < 190 GeV/c 2 . 

Measurements of the cross section for top quark pair production provide a test of the 
expected QCD production mechanism as well as of the standard model decay into a W- 
boson and a bottom quark, t — > Wb. Non-standard model production mechanisms could 
enhance the measured cross section, and non-standard model decays could suppress the 
measured value, which assumes a branching fraction of t — > Wb of 100%. 

In this paper we describe a measurement of the ti production cross section in pp collisions 
at = 1.96 TeV with the CDF II detector at the Fermilab Tevatron. The measurement 
assumes the standard model decay t — > Wb of the top quark, providing a final state from 
ti production that includes two W bosons and two bottom quarks. We select events where 
one of the W bosons decays to an electron or muon which has large momentum transverse 
to the beam direction (Pt) plus a neutrino. The neutrino is undetected and results in an 
imbalance in transverse momentum. This imbalance is labeled "missing E^" (I$t) because it 
is reconstructed based on the flow of energy in the calorimeter [2] . The other W boson in the 
event decays hadronically to a pair of quarks. The two quarks from the W boson, and the 
two b quarks from the top decays, hadronize and are observed as jets of charged and neutral 
particles. We take advantage of the semileptonic decay of 6-hadrons to muons to identify 
final-state jets that result from hadronization of the bottom quarks expected in the top decay. 
This "soft-lepton tagging" with muons, or SLT/i, is effective in reducing the background 
to the ti signal from W plus multijet production. This technique is complementary to 
measurements that take advantage of the long lifetime of 6-hadrons to identify jets from 
bottom quark hadronization through the presence of a decay vertex displaced from the 
primary interaction [3]. 

This measurement is an update of the measurement described in [4], which was made 
with approximately one tenth of the integrated luminosity used here. Full details of this 
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analysis are presented in [5]. In addition to the larger dataset, we report here on a new 
method for evaluating the background from "mistags" , i.e. those SLT/is that do not arise 
from semi-leptonic decays of heavy-flavor (HF) quarks. This is described in Section VI. 

II. THE CDF DETECTOR 

The CDF II detector is described in detail in [6]. We describe briefly here those elements 
of the detector that are central to this analysis. CDF II is a nearly azimuthally and forward- 
backward symmetric detector designed to study pp interactions at the Fermilab Tevatron. 
It consists of a magnetic spectrometer surrounded by calorimeters and muon chambers. An 
elevation view of the CDF II detector is shown in Fig. 1. 




FIG. 1: Elevation view of the CDF II detector. 

Charged particles are tracked inside a 1.4 T solenoidal magnetic field by an 8- layer silicon 
strip detector, covering radii from 1.5 cm to 28 cm, followed by the central outer tracker 
(COT), an open-cell drift chamber that provides up to 96 measurements of charged particle 
position over the radial region from 40 cm to 137 cm. The 96 COT measurements are 
arranged in 8 "superlayers" of 12 sense wires each, that alternate between axial and 2° 
stereo orientations. The silicon detector tracks charged particles with high efficiency for 
|77|<2.0, and the COT for |^|<1.0 [2]. 
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Surrounding the tracking system, and outside the magnet coil, are the electromagnetic 
and hadronic calorimeters, used to measure charged and neutral particle energies. The 
electromagnetic calorimeter is a lead-scintillator sandwich and the hadronic calorimeter is 
an iron-scintillator sandwich. Both calorimeters are segmented in azimuth and polar angle 
to provide directional information for the energy deposition. The segmentation varies with 
position on the detector and is 15° in azimuth by 0.1 units of rj in the central region (| 
rj |< 1.1). Segmentation in the plug region (1.1 <| rj |< 3.6) is 7.5° (| r\ |< 2.1) or 15° 
(| i] |> 2.1) in azimuth and ranges from 0.1 to 0.64 units of i] (corresponding to a nearly 
constant 2.7° change in polar angle). The electromagnetic calorimeters are instrumented 
with proportional and scintillating strip detectors that measure the transverse profile of 
electromagnetic showers at a depth corresponding to the shower maximum. 

Behind the central calorimeter are four layers of central muon drift chambers covering 
| rj |< 0.6 (CMU). The calorimeter provides approximately one meter of steel shielding. 
Behind an additional 60 cm of steel in the central region sit an additional four layers of muon 
drift chambers (CMP) arranged in a box-shaped layout around the central detector. Central 
muon extension (CMX) chambers, which are arrayed in a conical geometry, provide muon 
detection for the region 0.6 <| i] |< 1.0 with between four and six layers of drift chamber, 
depending on zenith angle. The CMX chambers covering from 225° to 315° in azimuth 
are known as the 'miniskirt' while those covering from 75° to 105° in azimuth are known 
as the 'keystone'. The remainder of the CMX chambers are referred to as the 'arches'. 
The muon chambers measure the coordinate of hits in the drift direction, x, via a drift 
time measurement and a calibrated drift velocity, and for CMU and CMX, the longitudinal 
coordinate, z. The longitudinal coordinate is measured in CMU by comparing the pulse 
heights, encoded in time-over-threshold, of pulses at opposite ends of the sense wire. In 
CMX, the conical geometry provides a small stereo angle from which the z coordinate of 
track segments can be measured. Reconstructed track segments in CMU and CMP have a 
maximum of 4 hits, and in CMX a maximum of 6 hits. 

III. DATA SAMPLE AND EVENT SELECTION 

This analysis is based on an integrated luminosity of 2034 ±120 pb" 1 [7] (1993 pb" 1 with 
the CMX detector operational) collected with the CDF II detector between March 2002 and 
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May 2007. 



A. Kinematic Selection 

The triggering and offline event selection used in this analysis are nearly identical to that 
used in the previous analysis described in [4]. For completeness we reproduce the basic 
trigger and selection criteria here and highlight the few differences. 

CDF II employs a three level trigger system, the first two consisting of special purpose 
hardware and the third consisting of a farm of commodity computers. Triggers for this anal- 
ysis are based on selecting high transverse momentum electrons and muons. The electron 
sample is triggered as follows: At the first trigger level, events are selected by requiring a 
track with Pp > 8 GeV/c matched to an electromagnetic calorimeter tower with Pt > 8 GeV 
and little energy in the hadronic calorimeter behind it. At the second trigger level, calorime- 
ter energy clusters are assembled, and the track found at the first level must be matched 
to an electromagnetic cluster with P T > 16 GeV. At the third level, offline reconstruction 
is performed and an electron candidate with P T > 18 GeV is required. The muon sample 
trigger begins at the first trigger level with a track with P T > 4 GeV/c matched to hits in 
the CMU and CMP chambers or a track with Pp > 8 GeV/c matched to hits in the CMX 
chambers. At the second level a track with Pp > 8 GeV/c is required in the event for all but 
the first few percent of the integrated luminosity, for which triggers at the first level were 
fed directly to the third level trigger. At the third trigger level a reconstructed track with 
P T > 18 GeV/c is required to be matched to the muon chamber hits. 

From the inclusive lepton dataset produced by the electron and muon triggers described 
above, we select offline an inclusive W plus jets candidate sample by requiring a recon- 
structed isolated electron with Pt > 20 GeV or muon with Pt > 20 GeV/c, ]$t > 30 GeV 
and at least 1 jet with Pt > 20 GeV and \rj\< 2.0. We define an isolation parameter, J, 
as the calorimeter energy in a cone of AR = a/ A?? 2 + A0 2 < 0.4 around the lepton (not 
including the lepton energy itself) divided by the E T (P T ) of the lepton. We select isolated 
electrons (muons) by requiring I < 0.1. Electrons and muons satisfying these criteria are 
called the "primary lepton". Jets are identified using a fixed-cone algorithm with a cone 
size of AP = 0.4 and are constrained to originate at the pp collision vertex. Their energies 
are corrected to account for detector response variations in r], drifts in calorimeter gain, 
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nonlinearity of calorimeter energy response, multiple interactions in an event and for energy 
loss in un-instrumented regions of the detector. These corrections bring the jet energies, on 
average, back to the sum P T of the particles in the jet cone, but not all the way back to the 
parton energy. This is slightly different from the previous analysis [4] where the correction 
was done only for response variations in rj, gain drifts and multiple interactions. The jet 
counting threshold in that analysis was > 15 GeV, which corresponds roughly to the 
E T > 20 GeV used here with the additional corrections. The missing transverse energy is 
corrected to account for the shifts in jet energies due to the jet corrections above, and the 
If/T threshold has been raised from 20 GeV in the previous analysis to 30 GeV here, consis- 
tent with the change in jet corrections. Z boson candidate events are rejected by removing 
events in which a second, same flavor, opposite sign isolated lepton, together with the pri- 
mary lepton, makes an invariant mass between 76 GeV/c 2 and 106 GeV/c 2 . The acceptance 
of these selection criteria for tt events is discussed in Section V below. 

The tt signal region consists of W candidate events with 3 or more jets, while the W+l 
and W+2 jet events provide a control sample with little signal contamination. 

The dataset selected above is dominated by QCD production of W bosons with multiple 
jets. As a first stage of background reduction, we define a total event transverse energy, H r ^ ) 
as the scalar sum of the electron E^ or muon Pp, J$t and jet for jets with Pp > 8 GeV and 
\rj\< 2.4. Figure 2 shows a comparison between the H T distributions for simulated tt and 
W+ jets events with at least 3 jets. For 3 or more jet events, we require H T > 200 GeV. This 
requirement rejects approximately 30% of the background while retaining approximately 
99% of the tt signal. No H T requirement is made for the control region of 1- and 2-jet 
events. 

B. Muon Tagging 

Even after the H T requirement is imposed, the expected signal to background ratio in 
W+ > 3 jet events is about 1:7. To further improve the signal to background ratio, events 
with one or more 6-jets are identified by searching inside jets for semileptonic decays of 
6-hadrons into muons. The algorithm for identifying such candidate 6-jets is known as "soft 
lepton tagging" or "SLT/i" and a jet with a candidate semileptonic b decay to a muon is a 
"tagged" jet. The SLT/i algorithm is described in detail in Reference [4]. We review here 
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FIG. 2: The distribution of for simulated tt and W+ jets events with at least three jets. 



only its basic features. 

Muon identification at CDF relies on the presence of a track segment ("stub") in the 
muon chambers, matched to a track in the central tracking system. The soft muon tagger 
is based on a x 2 function that uses all available information about the match between the 
extrapolated COT track and the muon stub to require that the deviations be consistent with 
the multiple Coulomb scattering expected for a muon traversing the CDF calorimeter. 

The algorithm begins by selecting "taggable" tracks. A track is declared taggable if it 
contains at least 3 axial and 2 stereo COT superlayers that have at least 5 hits each. To 
obtain some rejection for decays- in- flight (DIF), the impact parameter, do, of the track with 
respect to the beamline, is required to be less than 2 mm. The track is further required 
to originate within 60 cm of the center of the detector along the beam direction. Finally, 
the track must have a Pp above an approximate range-out threshold of 3.0 GeV/c and 
extrapolate to within a fiducial volume at the muon chambers that extends Sctms outside 
of the physical edges of the chambers, where <tms is the deviation expected from multiple 
Coulomb scattering at the track Pp. 

Matching between the extrapolated COT track and the muon stub is done using the 
following observables ("matching variables"): The extrapolated position along the muon 
chamber drift direction (x), the longitudinal coordinate along the chamber wires (z) when 
such information is available, and the extrapolated slope (4>l)- Tracks in the COT are paired 
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with stubs based on the best match in x, which must be less than 50 cm for a track-stub 
pair to become a muon candidate. We refer to the difference between the extrapolated and 
measured positions in x and z as Ax and Az, respectively, and between the extrapolated 
and measured slope as A<pi. The distributions of these variables over an ensemble of events 
are referred to as the "matching distributions" . 

Candidate muons are selected with the SLT/x algorithm by constructing a global x 2 quan- 
tity, L, based on a comparison of the measured matching variables with their expectations. 
The first step in constructing L is taking a sum, Q, of individual x 2 variables: 



where [ii and <7j are, respectively, the expected mean and width of the distribution of the 
matching variable JQ. The sum is taken over n selected variables as described below. We 
construct L, by normalizing Q according to 



where the variance, var(Q), is calculated using the full covariance matrix for the selected 
variables. The normalization is chosen to make L independent of the number of variables n. 

The selected variables are the full set of matching variables, Ax, Az, m the CMU, 
CMP and CMX with the following two exceptions: The CMP chambers do not provide 
a measurement of the longitudinal coordinate z, and matching in <f> L is not included for 
stubs in the muon chambers that have only three hits. Because of their significantly poorer 
resolution, track segments reconstructed only in the CMU or only in the CMP chambers 
with only three hits are rejected (if the SLT/z candidate has stubs in both CMU and CMP, 
then a stub with only three hits is allowed). These two exceptions are a new feature of the 
algorithm, since the previous publication, that reduce backgrounds from hadronic punch- 
through with a negligible effect on the efficiency Note that a muon that traverses both the 
CMU and the CMP chambers yields two sets of matching measurements in x and 4>l and 
one z matching measurement, and is referred to as a CMUP muon. All available matching 
variables are used in the calculation of L for a given muon candidate. 

As described in Reference [4], the expected means and widths in equation 1 are param- 
eterized as a function of the Pt of the muon using J ftp and W and Z bosons in the data. 




(1) 



L = 



(Q-n) 



(2) 



^var(Q) ' 
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We use the same parametrization described there. The efficiency has been remeasured, from 
the data, using the full dataset for this analysis. 

Using J/ip events only we measure the efficiency as a function of the quantity L defined in 
equation 2 (the efficiency measurement is described in detail in Section VB). The efficiency 
plateaus at a value of | L | < 3.5, and we therefore use this requirement to define an SLT/i 
tag. 

Beginning with the lU+jets candidate dataset, selected as described in Section III A 
above, we require that at least one jet in each event has an SLT/i tag. A jet is determined 
to have an SLT/i tag if a candidate muon with | L |< 3.5 is found within a cone of AR < 0.6 
centered on the jet axis. When the primary lepton is a muon, the event is rejected when the 
SLT/i has opposite charge to the primary muon and together with that muon has an invariant 
mass between 8 and 11 GeV/c 2 or between 70 and 110 GeV/c 2 . This rejects events in which 
an T or Z boson decays to a pair of muons, one of which becomes the primary lepton while 
the other ends up in a jet and is tagged by the SLT/i algorithm. Whether the primary is an 
electron or a muon, events where the invariant mass is less than 5 GeV/c 2 are also removed 
to prevent sequential double-semileptonic b — > c — > s decays (where the primary lepton and 
the SLT/x tag are from these semileptonic decays, rather than the primary lepton being from 
the decay of a W boson) from entering the sample, as well as events with a J/ip decay. We 
further reject events as candidate radiative Drell-Yan and Z bosons if the tagged jet has an 
electromagnetic energy fraction above 0.8 and only one track with P T > l.OGeV/c within a 
cone of AR = 0.4 about the jet axis. 

Three levels of selection are defined in this analysis. Events that pass the kinematic 
cuts and the dilepton and radiative- Z vetoes, but do not necessarily have an SLT/i-taggable 
track in them, comprise the "pretag" sample. Pretag events that have an SLT/i-taggable 
track (P T > 3 GeV/c, passing quality cuts, pointing to the muon chambers) within AR < 0.6 
of a jet with Ei > 20 GeV are called taggable events. Finally, the subset of SLT/i-taggable 
events that have at least one SLT/i-tagged jet are called tagged events. 

C. Selected Event Samples 

Table I shows the number of pretagged, taggable and tagged events in the electron and 
muon channels in this dataset as a function of jet multiplicity. 
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1 jet 2 jets 3 jets >4 jets >3 jets 
Electrons 

Pretag 79348 13068 1615 660 2275 

SLT> Taggable 43005 10479 1518 648 2166 

SLT> Tagged 519 224 85 64 149 

CMUP Muons 

Pretag 38165 6320 719 325 1044 

SLTTi Taggable 20162 4921 673 312 985 

SKl> Tagged 224 105 41 34 75 

CMX Muons 

Pretag 23503 3672 422 162 584 

SM> Taggable 12428 2864 396 160 556 

SLTjU Tagged 149 55 16 8 24 

Electrons+Muons 

Pretag 141016 23060 2756 1147 3903 

SLTfi Taggable 75595 18264 2587 1120 3707 

SLT/i Tagged 892 384 142 106 248 

TABLE I: Summary of event counts for 2 fb _1 of CDF Run II data for the event selection described 
in Sections III A and IIIB. 

IV. MONTE CARLO DATASETS 

The detector acceptance of tt events is modeled using PYTHIA v6.216 [8] and 
herwig v.6.510 [9]. This analysis uses the former for the final cross section estimate and 
the latter to estimate the systematics resulting in the modeling of tt production and decay. 
The PYTHIA event generator has been tuned using jet data to better model the effects of 
multiple interactions and remnants from the break-up of the proton and antiproton. The 
generators are used with the CTEQ5L parton distribution functions [10]. Decays of b- and 
c-hadrons are modeled using evtgen [11]. 

Events with a W boson produced in association with multiple jets are modeled using 
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alpgen v2.1 [13], with parton showering provided by pythia v6.326 and HF hadron decays 
handled by evtgen. alpgen calculates exact matrix elements at leading order for a large 
set of parton level processes in QCD and electroweak interactions. The showering in pythia 
may result in multiple alpgen samples covering the same phase space. These overlaps are 
removed using a jet-parton matching algorithm along with a jet-based heavy flavor overlap 
removal algorithm [12]. 

Estimates of backgrounds from diboson production (WW, WZ, and ZZ) and Drell- 
Yan/Z — > rr are derived using pythia. Drell-Yan to /i/i events are modeled using alpgen 
with pythia showering while single-top production is modeled with madevent [14], also 
with PYTHIA showering. 

The CDF II detector simulation reproduces the response of the detector to particles 
produced in pp collisions. The detector geometry database used in the simulation is the 
same as that used for reconstruction of the collision data. Details of the CDF II simulation 
can be found in [15]. 

V. EFFICIENCY FOR IDENTIFYING ti EVENTS 

The efficiency for identifying ti events in this analysis is factorized into the geometric 
times kinematic acceptance and the SLT/i tagging efficiency The acceptance is evaluated 
assuming a top mass of 175 GeV/c 2 and includes the branching fraction to leptons, which 
is assumed to have the SM value. The tagging efficiency is the efficiency for SLT//-tagging 
at least one jet in events that pass the geometric and kinematic selection. Each piece is 
described below. 

A. Geometric and Kinematic Acceptance 

The acceptance of ti events in this analysis is measured in pythia and then corrected, 
using measurements from the data, for effects that are not sufficiently well modeled in 
the simulation: the lepton trigger efficiencies, the fraction of the pp luminous region well- 
contained in the CDF detector (i.e. the z- vertex cut efficiency), and track reconstruction 
and lepton identification efficiencies. The efficiency of the z- vertex cut, | z |< 60 cm, is 
measured from minimum-bias triggered events to be (96.3 ± 0.2)%. The correction factor 
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for the difference between the track reconstruction efficiencies in data and simulation is 
1.014 ± 0.002. Events in the Monte Carlo are not required to pass any trigger, so the 
acceptance is multiplied by lepton trigger efficiency. The lepton trigger and identification 
efficiencies, measured using the unbiased leg of Z boson decays to electrons and muons, and 
the correction factors for each of the primary lepton types are shown in Tab. II. 



Quantity 


Electron 


CMUP Muon CMX Muon 


Trigger efficiency 


0.966±0.005 


0.917±0.005 


0.925±0.007 


Lepton ID efficiency (data) 
Lepton ID efficiency (MC) 


0.789±0.004 
0.806±0.001 


0.829±0.006 
0.896±0.001 


0.893±0.006 
0.916±0.002 


Lepton identification correction 0.978±0.005 


0.926±0.007 


0.975±0.007 



TABLE II: Summary of lepton trigger and identification efficiencies. 



The raw acceptance is defined as the number of pretag events divided by the total number 
of ti events in the pythia sample. The acceptance, after correcting for the differences 
between data and simulation, is shown in Tab. Ill as a function of the number of identified 
jets. 

1 jet 2 jets 3 jets > 4 jets > 3 jets 

Electron (%) 0.163±0.002 0.858±0.004 1.63±0.01 2.08±0.01 3.71±0.01 

CMUP Muon (%) 0.088±0.001 0.472±0.003 0.909±0.004 1.142±0.005 2.05±0.01 
CMX Muon (%) 0.042±0.001 0.219±0.002 0.414±0.003 0.532±0.003 0.946±0.004 
Combined (%) 0.292±0.002 1.544±0.005 2.946±0.008 3.743±0.009 6.69±0.01 

TABLE III: Acceptance for ti events as a function of jet multiplicity from pythia Monte Carlo 
sample, after data/MC corrections described in the text. In the combined acceptance we account 
for the fact that the CMX detector was not operating early in Run II. The uncertainties listed are 
statistical only. 
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B. Efficiency of the SLT> Algorithm 



The efficiency of the SLTfx algorithm is measured from the data using samples of J/ip and 
Z decays triggered on a single muon. The tagger is applied to the non-trigger muon (probe 
leg). If both legs pass the trigger, only one of them is used. To reduce background in 
the Z sample, the leg that is not used to measure efficiency is required to be isolated and 
to be consistent with being a minimum ionizing particle in the calorimeter. We correct 
for the remaining background using the invariant mass regions outside the Z mass window 
("sidebands"). 

The efficiency of the SLT/x is defined as: 

Number of tagged muons 
Number of taggable muon tracks with a stub 

The requirement in the denominator that the taggable muon track has a stub in the 
requisite muon chambers decouples the muon reconstruction efficiency, which is accounted 
for separately, from the efficiency of the tagger. Figures 3 and 4 show the efficiency for 
tagging muons with | L \ < 3.5 as a function of muon Pp from both J/ip and Z data. 
The decrease in efficiency with increasing Pt is due, primarily, to non-Gaussian tails in the 
resolution functions. These efficiency data are fit to functional forms [4] shown as the curves 
in the data. The dotted curves are those obtained by varying the fit parameters by ±lcx. 
Although the efficiency measurement is dominated by isolated muons, we do not expect that 
it will depend on the isolation of the muon because the muon chambers are well shielded 
from the inner detector. We have checked this assumption by measuring the efficiency as a 
function of the number of nearby tracks and found no dependence. 

The efficiency for SLT/x-tagging a ti event is measured for Monte Carlo events that pass 
the geometric and kinematic selection. We model the SLT/z tagging in these events by 
tagging muons from semileptonic HF decay with a probability given by the efficiencies in 
Fig.s 3 and 4. Events without a 'real' muon tag of this type can still be SLT/x-tagged through 
a mistag. Mistags in ti events are included by applying the "mistag matrix" described in 
the following section, and are included as part of the signal efficiency. The SLT/z tagging 
efficiency in ti events is given in Tab. IV. 
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FIG. 3: The SLT// efficiency for CMU/CMP as a function of Pt measured from J/iJj and Z data 
for | L |< 3.5. The solid line is the fit to the data and the dashed lines indicate the uncertainty on 
the fit. 



95- 




p T [GeV/c] 



FIG. 4: The SLT/i efficiency for CMX Arches (circles) and Miniskirt /Keystone (triangles) as a 
function of Pt measured from J/ip and Z data for | L |< 3.5. The solid curves are the fits to the 
data and the dashed lines indicate the uncertainties on the fits. 

VI. PREDICTING THE NUMBER OF TAGS FROM LIGHT-QUARK JETS 

As a prelude to the evaluation of the backgrounds to the ti signal we describe a new 
method, developed for this analysis, for predicting the number of SLT/i tags that come from 
light quark jets. We refer to these as "mistags", and they result from a combination of 
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1 jet 2 jets 3 jets > 4 jets > 3 jets 

Electron (%) 7.0±0.3 11.4±0.2 12.9±0.1 14.9±0.1 14.0±0.1 

CMUP Muon (%) 5.6±0.3 10.7±0.2 11.8±0.1 14.1±0.1 13.1±0.1 

CMX Muon (%) 6.7±0.5 11.2±0.3 12.3±0.2 14.2±0.2 13.4±0.2 

Average (%) 6.5±0.2 11.2±0.1 12.5±0.1 14.6±0.1 13.6±0.1 

TABLE IV: tt event tagging efficiency for SLT muons as a function of jet multiplicity from pythia 
Monte Carlo sample. The lepton category refers to the primary lepton. The average tagging 
efficiency is determined by weighting each channel by the acceptance and luminosity for each 
channel. The listed uncertainties are statistical only. 

hadronic punch-through of the calorimeter and muon steel, and hadronic decays-in-flight. 

To predict the number of mistags in our sample we use a track-based mistag probability 
that is a function of track P T and r\. We use reconstructed D* and A to identify a clean 
sample of pions, kaons and protons and measure the probability per taggable track, in 8 bins 
of Pt and 9 bins of 77, for each to satisfy the SLT/z | L |< 3.5 requirement. Details of the 
reconstruction technique, the measurement of the tagging probabilities and the assembly 
and testing of the two-dimensional (8 x 9) "mistag matrix" are described in what follows. 

A. Data samples 

To identify kaons and pions we reconstruct D* + — > D°tt + — > K~7t + tt + decays, and their 
charge conjugates. This dataset is collected using a two-track trigger that requires two 
oppositely charged tracks with Pt >2 GeV/c. The tracks are also required to have a scalar 
sum P T1 + P T2 > 5.5 GeV/c, an opening angle between them of 2° < |A0| < 90°, and 
originate from a displaced vertex. 

A sample of protons is obtained by reconstructing A — > pn~ decays. These events are 
collected using another two-track trigger similar to the one described above, but with an 
opening angle requirement of 20° < | A0| < 135° and the invariant mass of the track pair (as- 
sumed to be pions) required to be 4 GeV/c 2 < M(ir, n) < 7 GeV/c 2 . 
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B. Event Reconstruction 

We apply the following track quality criteria in the reconstruction of both D* [16] and 
A [17] decays: 

• the number of COT axial superlayers with > 5 hits is > 3; 

• the number of COT stereo superlayers with > 5 hits is > 2; 

• the track has \z \ < 60 cm. 

The D* reconstruction then proceeds through the examination of the mass difference 
Am = m(K-n-n) — m(Kix) with the following criteria: 

• the kaon must have opposite charge to each of the two pions; 

• \Azq\ < 5 cm between any two tracks; 

• the soft pion from the D* — > D°n decay must have Pp > 0.5 GeV/c; 

• the kaon and pion from the D° decay must each have Pr > 2 GeV/c; 

• the kaon and pion tracks must have impact parameter, \do\ < 0.2 cm; 

• \m(Kn) - m(D°)\ < 0.03 GeV/c 2 ; 

• At least one of the tracks (K or n) from the D° must be SLT/x taggable (including 
having P T > 3 GeV/c). 

As shown in Fig. 5, a clean D* signal is obtained for the right-sign Am distribution. 
The reconstruction of A decays requires the following criteria: 

• the pion and proton must have opposite charge; 

• \Az \ < 2 cm between the two tracks; 

• the x 2 of the vertex fit must be < 10; 

• the vertex must have L xy > 0.5 cm, where L xy is defined as the projection onto the 
net momentum direction, in the r — <fi plane, of the vector pointing from the primary 
to the secondary vertex.; 
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FIG. 5: The m(Kinr) — m(Kn) distribution for Z)* ± — > D Qr K^,D Q — > K^ir^ candidates in different 



SLT^-track-Px bins. The line in each plot represents the fit to the sideband regions. 
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FIG. 6: The m(pir) distribution for A — > pir candidates in different SLT/i-track-Px bins. The line 
in each plot represents the fit to the sideband regions. 

• the proton Pt is greater than the pion Pp; 

• the pion must have Pp > 0.4 GeV/c; 

• the proton must have |g?o| < 0.2 cm; 

• the proton must be SLT/z taggable (including having Pr > 3 GeV/c). 
Figure 6 shows the invariant mass distribution in the pir mass hypothesis. 
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Region 


Mass Window (MeV/c 2 ) 


D* Signal 


142.421 < Am < 148.421 


D* Sidebands 


139.6 < Am < 141 or 




152 < Am < 162.5 


A Signal 


1109.683 < m < 1121.683 


A Sidebands 


1090 < m < 1105.683 or 




1125.683 < m < 1170 



TABLE V: Mass windows used in determining the status of a D* or A candidate. 

We define a signal region for D* and A decays as well as sideband regions for each. We 
measure the sideband-subtracted tagging probability for K, n and p tracks using events in 
the signal, corrected for the enhanced probabilities in the sidebands. The sideband regions 
have a higher SLT/i per track tag probability because they are enriched in HF as a result of 
the two-track trigger described above. The signal and sideband regions are given in Tab. V 
and the sideband subtraction is done using the fits [5] shown in Fig.s 5 and 6. The tag 
probabilities before and after sideband subtraction are shown as a function of Pt in Fig- 7. 
We note that there are systematic uncertainties due to the choice of fit functions and in 
particular the quality of the fits in the sideband regions. These systematics, and all others 
associated with the construction of the mistag matrix, are evaluated by testing the predictive 
power of the matrix on a variety of independent data samples, as described in Section VIII B. 

The mistag matrix is designed to predict SLT/i tags that arise from both hadronic punch- 
through and decays-in-flight. When a pion or a kaon from a D* decays in flight, the track 
may be poorly reconstructed causing the reconstructed mass to fall outside of the signal 
region defined in Tab. V. We measure the size of this effect using D* decays in a Monte 
Carlo sample and make a correction. The correction factor is calculated in three bins in P T 
(limited by the sample size of the Monte Carlo) and shown in Tab. VI. Full details of the 
calculation of the correction factor are given in [5]. 
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FIG. 7: The measured (triangles) and sideband-subtracted (circles) tag probabilities as a function 
of track for (a) pions, (b) kaons, and (c) protons. The uncertainties shown are statistical only. 
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Frac. DIF After Reco. DIF Tag prob. Corr. Factor 


[GeV/c] 


[%1 \%\ 

[/OJ [/OJ 


[/OJ 




3-4 


u.^±u in u.uo u.zu in u.u^i 


A A R 4- O n 


i oc _|_ n no 


7T 4-6 


0.25 ± 0.04 0.16 ± 0.03 


60.1 ± 2.3 


1.16 ± 0.08 


> 6 


0.20 ± 0.04 0.16 ± 0.04 


75.6 ± 2.6 


1.09 ± 0.14 


3-4 


0.99 ± 0.09 0.77 ± 0.08 


10.2 ± 1.3 


1.04 ± 0.02 


K 4-6 


0.65 ± 0.06 0.51 ± 0.06 


10.8 ± 1.3 


1.02 ± 0.02 


> 6 


0.39 ± 0.05 0.23 ± 0.04 


18.2 ± 1.8 


1.05 ± 0.02 



TABLE VI: Relevant numbers, from Monte Carlo simulation, in the determination of the decay- 
in-flight correction for D* decays. These include the fraction of taggable tracks from D*s that 
decay- in-flight (DIF), the same fraction after all reconstruction requirements that fall inside the 
D* signal window (shown in Tab. V) and the probability for tagging a decay- in- flight. The last 
column gives the correction factor that is applied to the measured tag probability to account for 
the bias against decays- in- flight. 

C. The Mistag Matrix 

At this point we have SLT/x tag probabilities for tracks from 7r, K and p, corrected for 
backgrounds (sideband-subtracted) and for a bias against ir and K decays-in-flight. What 
remains is to assemble these separate SLT/i tag probabilities into a full mistag matrix that 
can be used to predict the number of tags in light-flavor jets in P^+jets events. 

To assemble the final mistag matrix, we take a weighted sum of the individual n, K and 
p matrices as follows: 

Mij = W n ■ M% + W K -M« + W P - Mfj. (4) 

where is the entry in the i th Pt and j th rj bin of the final matrix and M-J, and 
M[j are the corresponding entries in the ir, K and p matrices. The weights W n = 71.9%, 
Wk = 15.6% and W p = 12.5%, are taken from the taggable-track particle content of light- 
quark jets in ALPGEN W+jets Monte Carlo. Figure 8 shows the final tag probability for the 
eight Pt bins (integrated over rf) and the nine 77 bins (integrated over P T ). The features 
in the r\ distribution are due to the profile of absorber in front of, and the coverage of, the 
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FIG. 8: The mistag probability per track as a function of track Pt and detector rj. The histogram 
binning matches that of the matrix. The uncertainties shown are statistical only. The structure 
in the rj distribution is an artifact of the profile of the absorber and geometrical coverage of the 
muon system. 

muon system. The average tag probability per track in the matrix is (0.41 ± 0.01)% per 
track. 

The uncertainties associated with the probabilities in the mistag matrix include uncer- 
tainties from the sideband subtraction, the decay-in-flight correction and the weighting of tt, 
K and p probabilities, to name just a few. To evaluate the overall systematic uncertainty on 
the number of light-quark tags predicted by the matrix, one possibility would be to carefully 
evaluate the size of each of these uncertainties. However, there is no straightforward way 
to do this. Instead, as described in detail in Section VIII B, we directly test the predictive 
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power of the matrix, using event samples acquired with jet triggers, and use the control 
samples of W + 1 and 2 jet events to further validate the technique and establish that we 
have not underestimated the size of the systematic uncertainty. 



VII. BACKGROUND EVALUATION 

The dominant background contribution to the ti signal in this analysis comes from mistags 
in W+jets events. Another smaller, yet still significant background comes from W bosons 
produced in association with heavy flavor (Wbb, Wcc and Wc). The estimate of the mistag 
background is described in Section VII A, while the W^+HF background estimate is described 
in Section VII B. 

Other backgrounds that can produce a W boson and an SLT/i tag that are not accounted 
for by the mistag matrix include dibosons (WW, ZZ, WZ), Z — > r + r~, single top, QCD 
multijet backgrounds including bb, and residual Drell-Yan (/i/i) events not removed by the 
dimuon removal. QCD and Drell-Yan backgrounds are measured using the data, as described 
in detail in Sections VII C and VII D below. The remaining backgrounds are estimated from 
Monte Carlo as described in Section VII E. We treat QCD independently of the calculation 
of mistags in I^+jets events because events that enter our sample by mimicking the signature 
of a W boson can have a significantly larger tag rate than true W events. The enhanced tag 
rate arises due to the contribution of bb events to the QCD background and because of the 
correlation between the tag rate and measured I$t m events in which the ~E$t arises from 
jet mismeasurement or semileptonic HF decay rather than from a neutrino in a W boson 
decay. In order to avoid double counting we correct the estimate of tags in W^+jets events 
by (1 — Fqcd), where Fqqd is the QCD multijet fraction in the W^+jets candidate sample. 

A. Mistags 

The background due to mistags is evaluated using the track-based mistag matrix de- 
scribed in Section VI. To predict the number of events from W+jets with at least 1 mistag, 
we apply the mistag matrix to all pretag events according to: 
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where the sum runs over each event in the pretag sample, and the product is over each 
taggable track in the event. V(Pti, rji) is the probability from the mistag matrix for tagging 
the i th track with parameters P-pi and rji. Note that the sum over the events in equation 5 
includes any tt events that are in the pretag sample. We correct for the resulting overestimate 
of the background at the final stage of the cross section calculation (see Section IX). It also 
includes VE+HF events, diboson events, etc. Therefore mistags from these backgrounds are 
included here. Tags from muons resulting from the decay of HF hadrons or W or Z bosons 
in these backgrounds are calculated separately using Monte Carlo simulations, as described 
in Sections VII B and VII E. To avoid any double counting the Monte Carlo estimates of 
the contributions from these backgrounds do not include any mistags. 

A fraction, Fqcd, of the events in the signal region are QCD events for which the back- 
ground is estimated separately. Therefore, we correct the prediction of equation 5 according 
to 

N%*« = (1 - F QCD ) ■ NZt g - (6) 

The background estimate from the application of the mistag matrix is shown in Tab. VII. 
We list here both the raw prediction and that corrected by (1 — Fqcd)- The calculation of 
Fqcd is described in Section VII C I. 

1 jet 2 jet 3 jets > 4 jets > 3 jets 
i\S tag 641±32 238±12 55.0±2.8 32.7±1.6 87.5±4.4 
622±31 226±12 53.0±2.7 31.4±1.6 84.5±4.3 

TABLE VII: Summary of background estimate from mistags in IE+jets events. These numbers 
include a contribution from tt events in the TE+jets sample that is removed in the final cross section 
calculation, as described in Section IX 



B. lE+Heavy Flavor 

The evaluation of background tags from the semileptonic decays of HF quarks in Wbb, 
Wcc and Wc events is done using the alpgen Monte Carlo program. We determine the 
fraction of VE+jets events that contain heavy flavor at the pretag level, F HF , and the tagging 
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efficiency, 6hf, for these events and then normalize the total to the number of W+jets events 
seen in the data. The final prediction of the number of tags from H^+heavy-flavor events is: 

NffF — (1 — FqcD — Father) ' N pretag • F HF ■ €HF- (7) 

where Fqcd is the fraction of QCD events in the pretag sample and F ot h er is the fraction 
of other, non-H^+jets backgrounds. As with the mistag prediction, correction for tt in the 
pretag sample is done as part of the final cross section calculation. 

This procedure is used because the theory cross sections for the Wbb, Wcc and Wc 
processes have large uncertainties, whereas the uncertainties on the fraction of events with 
heavy-flavor jets are smaller. This procedure follows that used in [3]. 

1. Heavy-Flavor Fractions & Tagging Efficiency 

The HF fractions of events in the J¥+jets sample are determined by measuring the frac- 
tions in Monte Carlo and then scaling those fractions by a multiplicative factor of 1.15±0.35, 
determined by comparing measured HF fractions in inclusive jet data with those predicted 
by ALPGEN. 

The ALPGEN HF fractions, broken down according the number of b- or c-jets, are shown 
in Tab. VIII. In addition to the uncertainty on the HF-fraction scaling, an additional uncer- 
tainty on the alpgen fractions is determined by varying the ALPGEN generator parameters 
such as Q 2 and the quark masses. 

The Monte Carlo is also employed to determine the efficiency for tagging a muon from a 
semileptonic heavy-flavor decay in ly+heavy-flavor events. As with the tt tagging efficiency 
described in Section V, tags are assigned based on the SLT/z tagging efficiency measured in 
the data (Fig. 3 and 4). The results are shown in Tab. VIII. Note that we do not include here 
the additional efficiency that arises from mistags in real HF jets, because this is included in 
the mistag evaluation given in Tab. VII. 

Armed with these HF fractions and tagging efficiencies, the number of tagged events from 
Vy+HF is evaluated according to equation 7 above. The results are given in the last line of 
Tab. VIII. 
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1 jet 2 jets 


3 jets 


> 4 jets 


> 3 jets 


Category 2 b 


Fhf (%) 


0.9±0.3 


1.8±0.7 


2.8±1.1 


2.0±0.8 




7.8±0.2 


8.4±0.2 


8.5±0.3 


8.4±0.2 


Category 1 b 


Fhf (%) 


0.7±0.3 1.4±0.5 


2.6±1.0 


3.0±1.1 


2.7±1.0 


£hf 


3.54±0.05 4.30±0.06 


5.5±0.1 


5.8±0.2 


5.53±0.09 


Category 2 c 


Fhf (%) 


1.3±0.5 


2.8±1.1 


4.5±1.7 


3.1±1.2 


CHF 


3.1±0.1 


3.6±0.1 


3.5±0.2 


3.6±0.1 


Category 1 c 


Fhf (%) 


5.5±2.1 8.9±3.4 


11.0±4.1 


11.5±4.4 


11.1±4.2 




1.52±0.02 1.70±0.03 


2.04±0.07 


2.05±0.06 2.04±0.06 


Wbb+Wcc+Wc Background 


145±55 66.6±25.2 


15.3±5.8 


8.5±3.2 


23.0±8.7 



TABLE VIII: The heavy-flavor fractions, Fhf, tagging efficiencies, €hf, and VF+heavy-flavor 
background evaluated using alpgen Monte Carlo. The fractions are scaled by 1.15 as described in 
the text. The uncertainty on the heavy-flavor fractions includes that from the scaling factor and 
from variation of the ALPGEN parameters. 

C. QCD Background 

The background due to tags in QCD events that enter the signal sample is estimated 
by calculating the fraction of QCD events in the Vr+jets data and applying the standard 
mistag matrix times a multiplicative factor. The multiplicative factor is required because the 
tagging rate of QCD events that enter the pretag sample is higher than the corresponding 
tagging rate for PF+jets events. 
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1. The QCD Fraction 



The fraction of QCD events before SLT/i tagging is determined using the isolation, I 
(see Section III A), and I$ T of events with high-P T leptons and jets. Under the assumption 
that / and tf T are uncorrelated for QCD events, the number of QCD events in the ti signal 
region can be found by extrapolation from the non-signal regions 

NS CD = ^ (8) 
Where region D, the signal region, and regions A, B and C are defined according to 



Region A : $ T < 20 GeV 
Region B : $ T < 20 GeV 
Region C : $ T > 30 GeV 
Region D : $ T > 30 GeV 



/ > 0.2 
/ < 0.1 
/ > 0.2 
/ < 0.1. 



The event counts used in equation 8 are corrected for Monte Carlo predictions of the 
number of VT+jets and ti events in Regions A, B and C. The QCD fraction, Fq C d is then 
given by N® CD divided by the total number of events in Region D. The QCD fractions are 
given in Tab. IX. 

To evaluate the accuracy of the I$t- I prediction, two complementary regions in the plane 
are defined as: 

Region E : $ T < 20 GeV; O.K I < 0.2 
Region F : tf T > 30 GeV; 0.1 <I < 0.2. 

The different regions in the If?-! plane are shown in Fig. 9. Region F is outside the signal 
region and, once contamination from VF+jets and ti is removed, should have a QCD fraction, 
Fq CD , of approximately 1.0. Fq CD is given by: 

F _ N C -N E 

We use the difference of Fq CD from 1.0 to estimate a systematic uncertainty on the I^t vs- 
/ technique. The results are given in Tab. IX. Given the deviation from 1.0 in the > 3 jets 
data, we assign an 11% (120%) systematic uncertainty to Fqcd for electrons (muons). 
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FIG. 9: A diagram illustrating the regions denned in the I$t vs. I plane. 

1 jet 2 jets 3 jets > 4 jets > 3 jets 

Electron channel 

Fqcd 0.0423±0.0009 0.070±0.002 0.049±0.003 0.056±0.006 0.051±0.003 
Region F 0.95±0.04 0.97±0.06 0.84±0.10 1.06±0.24 0.89±0.09 

Muon channel 

Fqcd 0.0118±0.0004 0.020±0.001 0.013±0.004 0.007±0.004 0.011±0.003 
Region F 0.58±0.05 0.65±0.07 0.31±0.09 2.27±4.25 0.45±0.13 

TABLE IX: The fractions, Fqcd, of lepton-plus-jets events due to QCD multijet processes be- 
fore SLT/i tagging. The uncertainties on the Fqcd values are statistical only. Also shown is 
the measured QCD fraction in Region F, used to assign a systematic uncertainty on the Fqcd 
prediction. 



2. The Tag Rate of QCD Events 

The tag rate in QCD events that populate our signal region is enhanced relative to the 
rate predicted by the mistag matrix. There are two sources for this enhancement. First, 
much of the J^t in QCD events is due to mismeasurement of jet energies, which is correlated 
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FIG. 10: The ratio of observed to predicted tags as a function of I$t in data with a non-isolated 
primary lepton (/ > 0.2). 

with the tag rate (see Section VIII B). As seen in Fig. 10, the ratio of observed to predicted 
tags increases with I$t- Second, QCD includes bb and cc events in which the high-Pp lepton 
comes from the semileptonic decay of one of the b or c quarks; if the other also decays 
semileptonically, it may be tagged by the SLT/i. 

In the I$t-I plane, the region closest kinematically to the signal region (Region D) is the 
high-isolation and high-^y region, Region C. Region C has the same I$t requirement that 
Region D does, and likewise it requires a high-Pr lepton in the event. Therefore, the tag 
rate measured in Region C is a good representation of that of QCD events in Region D. 

The tag rate of QCD events is measured as an enhancement factor, k, times the PF+jets 
mistag probability. We calculate k as the ratio of observed to predicted (by the mistag 
matrix) SLTyU tags in Region C. The results are shown in Tab. X. 

3. The QCD Background Estimate 

Having determined Fq C d and k, the QCD background is given by: 

Nqcd = ■ k ■ Fqcd (10) 

where iV r ^ tag comes from equation 5. 
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1 jet 2 jets 3 jets > 4 jets > 3 jets 

Electron channel 
Fqcd 0.042±0.005 0.070±0.008 0.051 ± 0.006 

k 5.3±0.4 3.9±0.4 3.7 ±0.5 

Nqcd 82.6±11.9 38.1±6.1 6.3±1.2 3.6±0.7 9.9±1.9 

Muon channel 
0.0118±0.0004 0.0205±0.0009 0.011 ± 0.014 

2.9±0.4 3.4±0.4 3.0 ±0.5 

9.3±11.3 6.8±8.3 0.7±0.9 0.4±0.6 1.2±1.4 

Combined channels 
Total Nqcd 92±17 44.9±10.4 7.0±1.5 4.1±0.9 11.1±2.4 

TABLE X: Summary of the QCD background estimate. The uncertainties on the QCD frac- 
tions, Fqcd, and on the number of QCD events, Nqcd, are systematic and statistical combined. 
The Fqcd and k values in the third column apply to 3, > 4 and > 3 jets. 

D. Drell-Yan^ mi 

Drell-Yan— > \x\x events can survive the Z, T and J/ip vetoes if one muon leg fails the 
isolation requirement. We evaluate the number of residual Drell-Yan events that remain in 
the signal sample after the dimuon vetoes by measuring the number of Z — > fj,fj, events in the 
data, Nf , inside the Z mass window, where one leg of the Z is identified as an SLT/i (these 
events are normally removed from the signal sample). We then use an alpgen Z/Y — > fj,fx 
Monte Carlo sample to estimate the ratio, R out / m ^ of events outside the Z-mass window to 
events inside. To increase the statistical precision the ratio is measured with the ]$t and Ht 
and dilepton rejection cuts removed and without requiring that the SLT/z be inside a jet. 
We use the Monte Carlo to measure the ratio of efficiencies e out / e rn Q f these requirements. 



Fqcd 
k 

Nqcd 
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With these pieces, the number of residual Drell-Yan events, N DY is 

z a _ut e™'(# r ,ff r ,dilep,SLT-jet) 
DY taa n ' e in (tf T , H T , dilep, SLT-jet) ' 

The results of the Drell-Yan calculation are given in Tab. XI. 





1 jet 2 jets 


3 jets > 4 jets > 3 jets 


N t % (data) 


27 25 


3 3 


out 
R in 


0.325±0.003 0.315±0.005 


0.312±0.006 


e out (%) 


0.71±0.04 1.62±0.10 


2.85±0.24 


e in (%) 


1.03±0.05 3.10±0.12 


3.24±0.18 


out ^out 


0.223±0.018 0.165±0.013 


0.274±0.027 


Drell-Yan Total 


6.02±1.25 4.12±0.88 


0.82±0.44 0.00±0.19 0.82±0.48 



(11) 



TABLE XI: Drell-Yan background summary. Uncertainties are statistical only. The values in the 
third column apply to 3, > 4 and > 3 jets. 



E. Monte Carlo Driven Backgrounds 

Backgrounds from dibosons (WW, WZ, ZZ), Z — > t + t~, and single top are determined 
from Monte Carlo. For each of these backgrounds, the estimated number of tags is calculated 

as 

AT, 4- «./!*. (12) 

Here is the theoretical cross section. The acceptance, A i: and the SLT/i tagging efficiency, 
6j, are calculated from the Monte Carlo. As with the IF+HF evaluation in Section VII B, 
the efficiency includes only that due to tagging a muon from a semileptonic heavy-flavor 
decay. We do not include mistags in the efficiency evaluation because this is included as 
part of the background determined by the mistag matrix. The background evaluations are 
shown in Tab. XII. 

VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties in this analysis come from Monte Carlo modeling of the geo- 
metrical and kinematic acceptance, knowledge of the SLT/i tagging efficiency, the effect on 
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the acceptance of the uncertainty on the jet energy scale, uncertainties on the background 
predictions, and the uncertainty on the luminosity. The evaluation of the size of each of 
these uncertainties is described below. 

A. Systematic Uncertainties on Acceptance and Efficiency 

Monte Carlo modeling of geometrical and kinematic acceptance includes effects of parton 
distribution functions (PDFs), initial-state radiation (ISR), final-state radiation (FSR), and 
jet energy scale. These are estimated by comparing different choices for PDFs, varying ISR, 
FSR and the jet energy scale in the Monte Carlo and comparing the pythia generator with 

HERWIG. 

The PDF uncertainty is evaluated from 3 contributions. The first is obtained by varying 
the PDF according to the 20 CTEQ eigenvectors [10]to account for the uncertainty on the 
PDF fit. The second is the difference between the CTEQ5L PDF used for the acceptance 
measurement with that obtained using MRST98 [22] in the default configuration to account 
for the type of PDF fit used. The third is evaluated comparing the default MRST with two 
alternative choices of a s to get an estimate of the uncertainty due to the value of a s . The 
three contributions in quadrature yield an acceptance uncertainty of 0.9%. 

The uncertainty due to the limited knowledge of ISR is constrained by studies of radiation 
in Drell-Yan events in the data. We vary both ISR and FSR in the ti Monte Carlo within 
the allowed range and add the deviations in quadrature. The systematic uncertainty due to 
this effect is 0.8%. 

The uncertainty on the acceptance due to the uncertainty in the jet energy scale is 
measured by shifting the energies of the jets in ti Monte Carlo by ±la of the jet energy 
scale [18]. The resulting uncertainty on the acceptance is 4.1%. 

The effects of generator modeling of the ti kinematics are measured by comparing the 
acceptance from pythia and herwig. The result is a 2.4% uncertainty. 

As described in Section VA, a scale factor is applied to the ti Monte Carlo data set to 
correct for lepton ID efficiency differences between data and Monte Carlo. This scale factor 
has an associated uncertainty that yields a 2.9% uncertainty on the total ti acceptance. 

The systematic uncertainty on the SLT/z tagging efficiency in ti events is comprised of 
three parts. First, the uncertainty due to the Pt dependence of the SLT/z efficiency curves 
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which is evaluated by remeasuring the ti event tagging efficiency with the ±la curves shown 
as the dashed lines in Fig. 3 and 4. Next, the tagging efficiency measurement assumes that 
the efficiency for finding tracks in jets in the COT is properly modeled in the simulation. 
This assumption comes with a 5% uncertainty, which was evaluated by embedding MC 
tracks in data events. Finally, the statistical uncertainty on the SLT/z efficiency in tt events 
is absorbed as a systematic uncertainty. The three contributions combine to give a systematic 
uncertainty of 5.1%. 

Adding all these contributions in quadrature gives a total "Acceptance Modeling and 
Efficiency" systematic uncertainty of 7.7%. 

B. Systematic Uncertainty of the Mistag Prediction 

To measure the uncertainty on the predicted number of tags in light-flavor jets, we test 
the predictive power of the mistag matrix on large samples of events triggered on a single 
jet with an (uncorrected) E T threshold of 20, 50, 70 or 100 GeV. Care must be taken in 
several areas. Because of the E T threshold on one jet in the event, that 'trigger jet' will have 
a bias against SLT/i tags because particles that reach the muon chambers do not deposit 
all their energy in the calorimeter and therefore reduce the measured jet energy from its 
true value. Therefore, if no other jet is above the trigger threshold, we remove the trigger 
jet from the sample used to test the mistag prediction. If there is an additional jet above 
the trigger threshold, then all jets above threshold are used. The opposite effect occurs in 
jets that are measured well below trigger threshold. In a di-jet event triggered, for instance, 
with a 100 GeV threshold, a single recoil jet with energy well below 100 GeV is likely to be 
significantly mis-measured. Such jets have an enhanced rate of SLT/i tags relative to jets in 
VF+jets events because jet mis-measurement is correlated with the population of SLT// tags 
through detector cracks, hadronic punch-through of the calorimeter and real muon content. 
In addition to rejecting the trigger jet, we reject jets in di-jet events if the recoil jet falls 
below the trigger threshold. For events with higher jet multiplicities we use only tracks in 
jets that are separated from the trigger-jet axis by AR between 0.7 and 2.6. These various 
criteria have been chosen in order to provide, in the jet samples, a set of jets that are similar 
in terms of SLT/i tags to those found in W+jets events [5]. 

To increase the number of jets available for the study we use, in addition to the jet 
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triggered data, events triggered on a single photon candidate (7+jets) with a threshold of 
25, 50 or 70 GeV, and a EEr-triggered sample, triggered on a four-jet total energy of at 
least 100 GeV. Jets in the 7+jets events are selected in the same way as in the single jet 
triggered events. All jets above 20 GeV are used in the HE^ sample. 

Since the mistag matrix is designed to predict the number of tags from light-flavor jets, 
we must also suppress heavy-flavor jets in our sample. This is achieved by removing events 
in which any jet has an identified secondary vertex [3], or in which the mass of the tracks 
contained in a potential secondary vertex is greater than 0.3 GeV/c 2 , or in which any jet 
contains a track with an impact parameter significance (d /c<io) > 2. This is found [5] 
to provide sufficient suppression of heavy-flavor jets while leaving the remaining sample 
unbiased against decays-in-flight inside the jet. 

With the above jet selection, the systematic uncertainty is determined using the difference 
between the number of SLT/i tags predicted by the mistag matrix and those observed in the 
data. The results are shown, as a function of the E T of the jet in Tab. XIII. 

Finally, we use the E T spectrum of VT+jets events from the Monte Carlo to perform a 
weighted average over the deviations between predicted and observed tags given in Tab. XIII. 
The result is (Predicted SLT/i - Observed SLT/i)/Pred.=(0.1±4.4)%. We assign a systematic 
uncertainty of 5% on the prediction of the mistag matrix. 

C. Other Background Uncertainties 

1. W+ Heavy Flavor Uncertainties 

Three sources contribute to the uncertainty on the Wbb + Wcc+Wc background predic- 
tion: the choice of alpgen settings, the uncertainty associated with the scaling factor that 
takes the heavy-flavor fraction in alpgen to the data, and the uncertainty on the tagging ef- 
ficiency. The determination of the uncertainties on the alpgen settings and the scale factor 
are described in Section VII B. The alpgen settings contribute 23% to the W+heavy- 
flavor background uncertainty and the scale factor another 13%. The uncertainty on the 
heavy-flavor tagging efficiency is the same as that for the ti tagging efficiency described in 
Section VIII A. The correlation between the efficiency for the background determination 
and for the ti acceptance is taken into account. 
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2. QCD Background Uncertainties 

Uncertainties on the QCD background prediction are determined using the level of agree- 
ment between predicted and measured events in 'Region F', as described in Section VII C I. 
We assign a systematic uncertainty on the Fq C d measurement of 11% for electrons and 
120% for muons, given conservatively by the worst agreement of the Region F prediction in 
each case. We fold this in with the statistical uncertainty on the Fqqd determination, the 
uncertainty on the correction factor k, both given in Tab. X, and the 5% systematic uncer- 
tainty due to the application of the mistag matrix. The total QCD background uncertainty 
is 19% and 124% for electrons and muons, respectively. In the final determination of the 
QCD systematic, we add in quadrature the separate effects on the cross section of the QCD 
uncertainties for electrons and muons. The estimate of the QCD background is correlated 
with the estimates of the mistags and rU+heavy-flavor backgrounds (equations 6, 7, and 10). 
This is taken into account when determining the effect on the ti cross section. Together with 
a relatively small QCD fraction of the events, the result is a rather small effect on the cross 
section determination, despite the large uncertainty on the QCD fraction itself. 

3. Other Background Uncertainties 

The systematic uncertainty on the small Drell-Yan background is determined by the 
statistical uncertainty of the estimate. Uncertainties on the Monte Carlo background pre- 
dictions come from uncertainties in the cross sections for the various processes and from the 
event sizes of the Monte Carlo samples. This uncertainty is reflected in the uncertainties 
quoted in Section VII E. The combined uncertainty on the Drell-Yan and Monte-Carlo- 
derived backgrounds is 11%. 

D. Summary of Systematic Uncertainties 

The systematic uncertainties are summarized in Tab. XIV. An additional systematic due 
to the uncertainty on the luminosity determination (5.9% [7]), is treated separately. 
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IX. tt PRODUCTION CROSS SECTION 



Before calculating the cross section, the estimated number of background events is cor- 
rected for tt events in the pretag sample using a simple iterative procedure. This is required 
because we apply the mistag matrix to the events before tagging to estimate the mistag 
and QCD backgrounds and also use the pretag sample in the W+heavy-flavor background 
determination assuming no tt content. A summary of the number of observed events and the 
background predictions, both before and after the correction, as a function of the number 
of jets is given in Tab. XV. It is worth noting the excellent agreement between the expected 
and observed tagged events in the W+l and 2 jet samples, where the expectation is dom- 
inated by the mistag contribution, and the tt contribution is negligible. This is a further 
validation of the mistag matrix. 

The cross section is calculated as 

(13) 

A tt e tt J i-iOX 

where N obs is the number of events with > 3 jets in which at least one jet has an SLT^u 
tag, N bkg is the corrected background, A t j and e t i are the tt event acceptance and tagging 
efficiency, and J Ldt is the integrated luminosity. The acceptance and efficiency are discussed 
in Section V, and summarized for the signal region in Tab. XVI. We measure a total tt cross 
section of 

a(pp -> ttX) = 9.1 ± l.lt.H ± 0.6 pb, (14) 

where the first uncertainty is statistical, the second is systematic, and the third is from the 
luminosity. This cross section value uses acceptances and tagging efficiencies appropriate 
for a top mass of 175 GeV/c 2 . The acceptances and efficiencies, and therefore the calculated 
cross section, change slightly for other assumed top masses. The calculated cross section 
is 3% higher assuming a top mass of 170 GeV/c 2 , and 4% lower assuming a top mass of 
180 GeV/c 2 . As a check we also measure the cross section separately for events in which the 
primary lepton is an electron and in which it is a muon. We measure 9.5 ± 1.2 pb when the 
primary lepton is an electron and 8.5 ± 1.2 pb when it is a muon. The uncertainties in both 
cases are statistical only. 

Figure 11 shows, in bins of the number of jets in VT+jets candidates, the expected number 
of tagged background and tt (normalized to the measured cross section) events together with 
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the number of observed SLT/x tags. 



• Data 
I |tt(c, r 9.1 pb) 

I Mistags 
HI Wbb.Wcc.Wc 
I I non-W 
S Other Bckg 

Background+tt 

Background+tt errors 

Require H T >200 GeV for > 3 jets 

* 



1 jet 2 jets 3 jets > 4 jets 

Jet Multiplicity 

FIG. 11: The expected background and observed tags in W+l, 2, 3, and 4-or-more jets events. 
The expected ti contribution is normalized to the measured cross section. 

In Fig. 12 through 14 we examine a few kinematic features of the tagged events. In 
each case the data are compared to the expected backgrounds plus tt, normalized to the 
measured cross section. The agreement between data and expectation is good. The only 
slight exceptions are a few bins at low in the W+ > 3 jet events in Fig. 12, where 
the number of observed tags exceeds somewhat the expectation. This is consistent with 
the excess seen in the low E^ jet data in Tab. XIII, which is folded into the systematic 
uncertainty on the measurement. 

X. CONCLUSIONS 

Using 2 fb -1 of integrated luminosity collected by the CDF II detector, we have measured 
the total cross section for tt production in pp collisions with a center-of-mass energy, y/s = 
1.96 TeV. The measurement begins by selecting a dataset of VF+jets candidates. We 
separate signal from background by identifying candidate semileptonic decays of 6-hadrons 
into muons. This technique was first published in Reference [4]. This measurement is an 
update that uses ten times the amount of data of the previous measurement and a new 
technique for evaluating the dominant background (see Section VI) of misidentifying a jet 



1000 



CO 

-4— ■ 






800- 


> 




(D 




JD 




Q) 
CT> 


600 


D) 




CO 




•*—> 




H — 




O 


400- 






n 




E 




z 


200- 



44 



> 300 
O 

« 250 
w 

200 

CD 

§150 
° 100 

CD 
_Q 

E 50 

3 







W + 1 ,2 jets 



• Data 
Mistags 
I I Wbb.WcC.Wc 
QCD 

Other backgrounds 
tt scaled to a . = 9.1 pb 
Background+tt uncertainties 




20 40 60 80 100 120 140 160 180 200 

Tagged Jet E T [GeV] 

W + > 3 jets 



Data 

tt scaled to o . = 9.1 pb 
Mistags 
I I Wbb,WcC,Wc 
QCD 

Other backgrounds 
Background+tt uncertainties 




20 40 60 80 100 120 140 160 180 200 

Tagged Jet E T [GeV] 



FIG. 12: Comparison of the jet Et distributions for tagged jets and for expectations from mistags, 
VF+heavy-flavor, QCD and events. The upper plot is for W+l- and 2-jet events and the lower 
plot for VF+3-or-more-jets events. 



from a light-flavor quark as one containing a 6-hadron. The measured tt cross section is 

a(pp -> ttX) = 9.1 ± l-l±5;g ± 0.6 pb, (15) 

consistent with the expectation of 6.7^g'g pb for standard model production and decay of 
top quark pairs with a mass of 175 GeV/c 2 . The measurement agrees well with other CDF 
measurements of the tt production cross section [23]. Assuming the cross section increases 
0.2 pb for every 1 GeV/c 2 decrease in the top mass, then at the world average top mass of 
172.4 GeV/c 2 the theoretical cross section is approximately 7.2 pb. Using a linear fit to the 
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FIG. 13: -Pp of the SLT// tags compared with expectations from backgrounds and tt. The upper 
plot is for W+l- and 2-jet events and the lower plot for W+ > 3-jet events. 



mass dependence, the measured cross section was estimated at the world average top mass 
and is found to be 8.9 ±1.6 pb. The kinematic distributions of the tagged sample are also 
consistent with standard model expectations. The observed number of tags in W + 1- and 
2-jet events is in excellent agreement with expectations from background, indicating that 
the backgrounds are well understood. 
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FIG. 14: The distribution of Pt relative to the jet axis (Pf rel ) for tags in data, compared with 
expectations from backgrounds plus tt. The upper plot is for 1- and 2-jet events and the lower plot 
for > 3-jet events. 
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1 jet 


2 jet 


3 jets 


> 4 jets 


> 3 jets 


WW 


^theory 






12.4±1.2 pb 






A{%) 


2.44±0.01 


2.62±0.01 


0.403±0.004 


0.121±0.002 


0.524±0.005 


e(%) 


0.49±0.03 


0.76±0.04 


0.88±0.09 


1.56±0.23 


1.06±0.09 


N ww 


2.986±0.299 5.001±0.394 


0.892±0.190 


0.475±0.118 


1.395±0.228 


WZ 


^theory 






3.96±0.40 pb 






A{%) 


1.085±0.007 1.317±0.007 


0.233±0.003 


0.070±0.002 


0.302±0.004 


e(%) 


0.85±0.06 


1.72±0.07 


1.46±0.16 


2.69±0.39 


1.77±0.15 


N W z 


0.740±0.075 1.821±0.128 


0.274±0.044 


0.151±0.036 


0.432±0.058 


ZZ 


^theory 






3.4±0.3 pb 






A{%) 


0.104±0.002 0.097±0.002 


0.060±0.001 


0.012±0.001 


0.042±0.001 


e(%) 


1.0±0.2 


2.4±0.3 


2.3±0.5 


1.6±0.5 


2.1±0.4 


iVzz 


0.07±0.02 


0.16±0.05 


0.05±0.02 


0.013±0.006 


0.06±0.02 






Drell - 


- Yan — > rr 






^theory 






333±4.2 pb 






A{%) 


0.112±0.001 0.054±0.001 0.0058±0.0004 0.0014±0.0002 0.0073±0.0004 


e(%) 


0.4±0.1 


0.4±0.1 


1.6±0.7 


1.3±0.5 


1.3±0.5 


^Drell-Yan->TT 


2.65±0.57 


1.54±0.43 


0.65±0.28 


0.13±0.05 


0.65±0.27 


s-channel Single Top 


^theory 






0.88±0.11 pb 






A{%) 


1.12±0.01 


2.66±0.01 


0.717±0.005 


0.203±0.003 


0.920±0.006 


e(%) 


5.0±0.1 


9.7±0.1 


10.2±0.2 


11.1±0.4 


10.4±0.2 


-^Ys— chan 


1.00±0.11 


4.61±0.46 


1.31±0.15 


0.40±0.05 


1.71±0.19 


t-channel Single Top 


^theory 






1.98±0.08 pb 






A{%) 


1.91±0.01 


2.10±0.01 


0.345±0.003 


0.057±0.001 


0.402±0.004 


e(%) 


4.38±0.09 


5.19±0.09 


5.96±0.24 


7.37±0.64 


6.16±0.22 


^Vt-chan 


3.36±0.37 


4.39±0.47 


0.83±0.11 


0.17±0.03 


1.00±0.13 



TABLE XII: Summary of Monte Carlo derived backgrounds. The theoretical cross sec- 
tion c M Ql [90l T91 1 q rp inplncnT-p TTi ci C'ccnr\i- q n r>& A inrOnrltiG f Vita Kr^rmViincr fr^c^f \r\r\ f n cnrcm^ a 



Jet E^ rr - 


\GeV] Observed 


Prodi rtoH 




20-30 


1 892 


i R41 + 29 


1 C O _|_ o o 






iuyo '-to 




45-65 


701 


768 ± 46 


8.7 ± 6.4 


65-90 


464 


462 ± 46 


-0.5 ± 11.0 


>90 


466 


466 ± 76 


-0.1 ± 16.9 


>20 


5084 


5029 ± 219 


-1.1 ± 4.6 



TABLE XIII: Checks of the mistag matrix in different jet E T bins. A =(Pred.-Obs.)/Pred. These 
values of A are weighted using the M^+3-or-more jets distribution to determine a systematic un- 
certainty on the mistag prediction. 



Source 


Fractional Sys. Uncert. 


(%) Aa ti (%) 


Acceptance Modeling and 
SLT Tagging Efficiency 


7.7 


+8.3 
-7.5 


Mistag Matrix Prediction 


5 


3.6 


Wbb + Wcc+Wc Prediction 


38 


5.3 


QCD Prediction 


19 (e) 124 (n) 


1.1 


Drell-Yan and other MC backgrounds 


11 


0.4 


Total Systematic Uncertainty 




+10.5 
-10.0 



TABLE XIV: Summary of systematic uncertainties. 
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H T > GeV H T > 200 GeV 



Background 


1 jet 


2 jet 


3 jets 


> 4 jets 


> 3 jets 


Taggable events 


75595 


18264 


2587 


1120 


3707 


Mistags 


622±31 


226±12 


53.0±2.7 


31.4±1.6 


84.5±4.3 


Wbb+Wcc+Wc 


145±55 


66.6±25.2 


15.3±5.8 


8.5±3.2 


23.0±8.7 


QCD multijet 


91.9±16.5 


44.9±10.4 


7.0±1.5 


4.1±0.9 


11.1±2.4 


ww+wz+zz 


3.80±0.44 


6.98±0.66 


1.21±0.23 


0.64±0.14 


1.88±0.30 


Drell-Yan^ t+t~ 


2.65±0.57 


1.54±0.43 


0.65±0.28 


0.13±0.05 


0.65±0.27 


Drcll-Yan^ 


6.02±1.25 


4.12±0.88 


0.82±0.44 


0.00±0.19 


0.82±0.48 


Single top 


4.36±0.39 


9.00±0.66 


2.14±0.18 


0.57±0.06 


2.71±0.23 


Total Background 


876.5±53.6 


359.0±24.0 


80.2±5.4 


45.3±3.0 


124.6±8.2 


Corrected Background 






79.5±5.3 


79.5±5.3 


ti Expectation (a =6.70) 


2.60±0.33 


23.5±1.8 


50.1±3.6 


74.2±6.5 


124.3±9.1 


Total Background + ti 


879.1±53.6 


382.5±24.1 


203.9±10.6 


203.9±10.6 


Tagged events 


892 


384 


142 


106 


248 



TABLE XV: Number of tagged events and the background summary. The uncertainty on the 
total background is not a simple sum in quadrature of the individual backgrounds because of the 
correlation between the mistag, VF+heavy-flavor, and QCD background predictions. 
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Electrons 


CMUP Muons CMX Muons 


Acc. no Tag (%) 


3.71 ±0.01 ±0.21 


2.05 ±0.01 ±0.14 0.946 ± 0.004 ± 0.050 


Event Tagging Eff. (%) 


14.02 ± 0.08 ± 0.72 


13.07 ±0.10 ±0.67 13.38 ±0.16 ±0.68 


Acc. with Tag (%) 


0.520 ± 0.003 ± 0.039 


0.268 ± 0.002 ± 0.022 0.127 ± 0.002 ± 0.009 


Luminosity (pb _1 ) 


2033.6 ± 119.6 


2033.6 ± 119.6 1992.5 ±117.2 


Denominator (pb _1 ) 


10.58 ±0.07 ±0.80 ±0.62 


7.97 ± 0.06 ± 0.49 ± 0.47 


Total denominator (pb _1 ) 


18.56 ± 0.09(stat.)±0.94(sys.)±1.09(lum.) 



TABLE XVI: Summary of components of the denominator for the cross section calculation. The 
ti acceptance and tagging efficiency for 3-or- more-jets events is determined using pythia Monte 



Carlo. 



53 



